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Abstract 

Soft-mode turbulence (SMT) is the spatiotemporal chaos observed in homeotropically aligned nematic 
liquid crystals, where non-thermal fluctuations are induced by nonlinear coupling between the Nambu- 
Goldstone and convective modes. The net and modal relaxations of the disorder pattern dynamics in SMT 
have been studied to construct the statistical physics of nonlinear nonequilibrium systems. The net re- 
laxation dynamics is well-described by a compressed exponential function and the modal one satisfies a 
dual structure, dynamic crossover accompanied by a breaking of time-reversal invariance. Because the 
net relaxation is described by a weighted mean of the modal ones with respect to the wave number, the 
compressed-exponential behavior emerges as a superposition of the dual structure. Here, we present exper- 
imental results of the power spectra to discuss the compressed-exponential behavior and the dual structure 
from a viewpoint of the harmonic analysis. We also derive a relationship of the power spectra from the 
evolution equation of the modal autocorrelation function. The formula will be helpful to study non-thermal 
fluctuations in experiments such as the scattering methods. 

PACS numbers: 05.45.-a, 61.30.-v, 47.54.De, 05.40.-a 
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I. INTRODUCTION 



Nonlinear phenomena are ubiquitous in nature and have various aspects. An example is chaotic 
(synonymously, turbulent) dynamics, in which few degrees of freedom (d.f.) generate unpre- 
dictable behavior (e.g., Refs. In spatially extended systems where the number of d.f. 
increases with expanding the system size, weak nonlinearity can lead to spatial and temporal dis- 
order. The phenomenon triggered by weak nonlinearity is called spatiotemporal chaos. There are 



several theoretical 
(KS) equation fa 



descriptions (Ref. yfl and references therein) such as the Kuramoto-Sivashinsky 
, the coupled map lattices [|5D, the Swift-Hohenbarg equation [|61], the Niko- 
laevskii equation [70, and the complex Ginzburg-Landau equation [|8|]. Nevertheless, systematical 
understanding for the chaotic or turbulent disorder is still developing. 

An experimentally obtained spatiotemporal chaos in homeotropically aligned nematic liquid 
crystals is soft-mode turbulence (SMT) [|9l-lllll. induced by nonlinear coupling between the local 
convective mode and global Nambu-Goldstone one (see Appendix for details). SMT has a char- 
acteristic spatial structure. Because a patchwork-like structure appears on the x-y plane when the 



angle of the local convective mode is color-coded, the structure is called patch structure [112 
The characteristic size ^ of the single patch domain is several of convective rolls and decreases as 
the distance from the onset of SMT as ^ ~ e^^/^, where e = (V/Vcf' — 1 is the normalized voltage 
used as a control parameter for SMT, V denotes the magnitude of the applied AC voltage, and 
the threshold voltage of the electroconvection. 

We have investigated SMT dynamics by measuring net temporal autocorrelation function Q{t) 
of turbulence-like dynamics in a steady state. It had long been considered that the simple exponen- 
tial described the relaxation dynamics. However, we revealed that the relaxation deviates from the 
simple exponential at the vicinity of the SMT's ons et 11 1511 : instead, it is well-fitted by the so-called 
Kohlrausch-Williams-Watts (KWW) function 1 1^ 

Q{t) = a exp 



(1) 



employed to explain the relaxation dynamics in glass-forming liquids (GFLs). It is considered that 
the cooperative rearranging in GFLs leads to the KWW-type relaxation I18]; indeed, spatially and 
temporally fluctuating domains have been experimentally observed near the glass transition point. 
The phenomenon is called the dynamical heterogeneity [Il9l-l22n. The characteristic length of the 
dynamical heterogeneity increases as the system approaches the glass transition point. Although 
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there are many studies for GFLs showing (3 < 1 (i.e., stretched exponential), the compressed- 



exponential behavior (i.e., (3 > 1) has been observed in some GFLs (e.g., Refs. |123 



m). In 



particular, Caronna et al. reported that the KWW exponent /3 increases toward 2 with approaching 
the glass transition point [26]. In SMT, the exponent is unity (i.e., simple exponential) at a large 
e and increases toward 2 (i.e., Gaussian) with decreasing e. The non-exponential behavior of 
the SMT's net relaxation originates from the patch structure. We have thus proposed a similarity 
between dynamics of SMT and GFL from a viewpoint of the spatial structure nl5\\ . 

The temporal correlations of each wave number, i.e., modal autocorrelation functions Ukir), 
are also suitable for studying chaotic or turbulent dynamics. Mori and Okamura have numerically 
and analytically discovered a dual structure in the KS equation the relaxation dynamics is 
dominated by the deterministic and stochastic orbits in the short-time and long-time region, re- 
spectively. We experimentally found the dual structure of the modal autocorrelation in SMT and 
have theoretically revealed the mechanisms of the dual structure ll29ll . It has been reported that 



the spatiotemporal disorder in SMT generates non-thermal fluctuations [|12l . 130(1 by which a non- 
Markovian (i.e., memory) effect is expected to play a significant role in the relaxation dynamics. 
We thus derived an evolution equation for the modal autocorrelation function in the projection- 
operator method for chaos and turbulence 113 ill and have specified the memory effect due to non- 
thermal fluctuations. The numerically obtained memory function from the experimental results 
has suggested that the non-thermal fluctuations are separated into Markov and non-Markov parts, 
where the latter part was named turbulent fluctuations. The relaxation dynamics is consequently 
separated into three relaxation stages: bare-friction stage, early stage, and late stage. The bare- 
friction stage disappears when the dissipation by the turbulent transport has much influence com- 
pared to that by the molecular one; subsequently, the modal relaxation in SMT is represented 
as 



' 1 - {r/i 



.(a) 



(V) 

{t <^tI ': early stage) 



Uk{r) oc < 



(2) 



exp 



-T / r. 



{rjf^ <^ r: late stage) 



where ' denotes the characteristic timescale of the memory due to the turbulent fluctuations. 
Furthermore, we have proved that the memory effect due to the turbulent fluctuations originates 
from the patch structure [|29D. 

In this proceedings paper, we uniformly argue the net and modal relaxation dynamics of SMT 
and discuss the compressed exponential behavior and the dual structure from a viewpoint of the 



harmonic analysis. This is useful for experiments such as the scattering methods. The contents are 
organized as follows. In Sec. |Ill the net and modal autocorrelation functions are defined and the 
memory function is reviewed according to Ref. n2% . In Sec. UlIl we show our experimental results 
of the power spectra and present the analytical relationship. We conclude this proceedings paper 
in Sec. |IVl In Appendix, the detail of SMT and the numerical algorithm for solving the memory 
function are summarized. 



n. RELAXATION DYNAMICS 



A. Net and modal correlation functions 



We observed the pattern dynamics of SMT through the transmitted light intensity /(r, t) with 
f = {x,y). Let Q{r,T) denote the autocorrelation function of the transmitted light intensity, 
defined as 

• ' (A/(r-,()2) • ^ ' 

where AI(f, t) = I{f, t) — {I{r, t)) denotes the fluctuation of the transmitted light intensity and 
the angle brackets indicate the long-time average in steady state; 
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T 



{fit + r)g{t)) = lim — / dtf{t + T)g{t). (4) 

T-i-OO 21 J _rp 

The net autocorrelation function Q{t) can be defined as the spatial average of Q{f, r); 

Q(r) = -^ / dfg(f,r), (5) 

Vd JVi 

where Vd denotes the generalized volume of rf-dimensional system and the integral range is whole 
system. The net relaxation in SMT is well described by the KWW function and the KWW expo- 
nent /3 is 1 at a large e and approaches to 2 with decreasing e [{l5^ . 

In addition, we have focused on the modal element of the fluctuation A/(r, t); 

uj:{t)= [ dfA/(r,t)e'^"^ (6) 

where i = Because SMT is isotropic [|9|], it is sufficient to study Uk{t), where k denotes the 

radial wave number; k = \k\. The modal autocorrelation function is defined as 

Mr) = <^i(^±lMM. (7) 



Here, the asterisk symbol denotes taking the complex conjugate and Pk = {\uk{t)\'^) the spatial 
power spectrum. The modal correlation function ?7fc(r) is a real number due to the translational 
symmetry and isotropy. In the rest of this paper, the wave number k is normalized by Aq as 
k = fcAo/27r, where Aq = Ao(£:) denotes the length of a pair of electroconvections. 

One can derive the analytical relationship between the net and modal autocorrelation functions. 
The autocorrelation of reduces to 



J fy^ dfdP (A/(r, t + r)A/(P, t)) e^^^-^^'^^. 



Integrating the both sides with respect to A;, we obtain 



For d = 2, the net autocorrelation function connects to the modal one as 



(8) 



Ah 



dA; kPkUkir) 

Q{r) = '-^-^ . (10) 

dA; kPk 



Therefore, in SMT, the compressed-exponential behavior appeared in the net relaxation is inter- 
preted as a superposition of the dual structure in the modal ones. 



B. Memory function 

We have proposed an approach by the memory function derived in the projection- operator 
formalism |29|] . The evolution equation for the modal autocorrelation function is represented as 



dUk{r) 
dr 



drriir-r')U,iT'), (11) 



where r^(r) denotes the memory function that results from the non-thermal fluctuations. The 
numerically-obtained memory function r^(r) has a sharp peak at r = 0, implying that the non- 



thermal fluctuations in SMT can be divided into Markov and non-Markov parts 



290,1.6., 



n(r) = 27r5(r) + r,(r). (12) 



where 7^°^ denotes the bare friction due to the Markov part of the non-thermal fluctuations and 
Ffc (r) the memory function caused from the non-Markov part. The evolution equation (fTT]) reduces 
to 

= -it^Ur) - /'drn(r - r')f/.(r'). (13) 

The non-Markov part of the non-thermal fluctuations is called the turbulent fluctuations. The 
algorithm to numerically solve the memory function is summarized in Appendix. 



The macroscopic friction 7^^^^ due to the turbulent fluctuations is represented as lUJ 



(r) 

Tk 



/ drrfc(r). (14) 
Jo 

Moreover, one can define a characteristic time scale of the memory function as 

' ~Jo ^^r,(o)-r,(o)' ^'^^ 

within which the memory effect due to the turbulent fluctuations is alive; in other words, for 

(r) 

T^, <ti T (late stage), the memory function due to the turbulent fluctuations is approximately 
represented by 

rfe(r) ^ 27^ (5(r). (16) 

fr) IT) 
In contrast, the non-Markovian regime r <^ ' can be analytically divided into r <^ Ik^k 

(bare-friction stage) and 'yk'^jf^ ^ r ^ rj^^ (early stage) |29], where 7^ denotes the ratio of the 

bare friction to the turbulent friction, 

7. = 7f /iP . (17) 

In the case that 7^ ^ 1, the bare-friction stage shrinks; hence, the early and late stages are observed 
as the dual structure in SMT. 



ni. HARMONIC ANALYSIS 

We first present the power spectrum of the fluctuation of /(r, t) in Fig. [H where Q^j was 
calculated from Q{t) as 

POO 

Qu; = 2 dTQ{T)cosujT. (18) 
Jo 

In the case that the relaxation satisfies the simple exponential function (i.e., /3 = 1 in the KWW 
function), the power spectrum is analytically derived as 



n — I — I — I I I 



n 1 1 — I — I — r—TT 




_] I I I I I 
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FIG. 1. Log-log plots of the power spectra eQuj (black line) for several control parameters; e = 0.050 
(top left), 0.10 (top right), 0.20 (bottom left), and 0.30 (bottom right). A gray solid line indicates the 
Lorenzian ( [20b and a gray dotted line the Gaussian (|2TI ). where A is set as 0.48 s^^. 

where Tq = TQ{e) denotes the characteristic correlation time of the relaxation. According to 
Ref. lllin . the above relationship reduces to 

2A 



A2 + n2 



(20) 



with A ^ = tqE and Q, = u/e. As discussed in Ref. llllh . the right-hand side of Eq. (|20|) does 
not depend on the control parameter e because tq is inversely proportional to e ||9|]. On the other 
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FIG. 2. Plot of the power spectrum Uf^^^ at e = 0.10 for several wave numbers; k = 0.76, 1.0, and 2.0 
from bottom to top. The inner plot is a closeup in an intermediate regime. 

hand, when the relaxation is the Gaussian (i.e., /3 = 2 in the KWW function), the power spectrum 
is analytically represented as 



eQui = exp 



4A2 



(21) 



To the best of our knowledge, no analytical functions for the Fourier transform of the KWW 
function with 1 < /3 < 2 have been proposed yet. As shown in the top panels of Fig. [H eQ^ for 
a small e, where (3 > 1, is in domains bounded by Eqs. (l20l) and (|2T1) . Morishita has proposed 
that the CEF-type decay appears as coexistence of ballistic and diffusive motion J32I. This plot 
style, eQuj versus Q, helps us to intuitively understand the KWW exponent 1 < (3 < 2. Further, as 
shown in the bottom panels of Fig. [T} Eq. (|20l) well describes the power spectrum at a higher e. It 
is evident that (3 converges to unity with increasing e. 

Next, we show results of the power spectrum U^^ui of the modal elements in Fig. |2l where Uk,ui 
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FIG. 3. Plot of the power spectrum at e = 0.10 for several wave numbers; k = 0.76, 1.0, and 2.0 
from bottom to top. 

was calculated by 

/•oo 

Uk,u^ = 2 dTUkiT)cOSUJT. (22) 



The harmonic analysis emphasizes on the short-time dynamics. Indeed, the results support appear- 
ance of the early and bare-friction stages in SMT dynamics: As shown in the inner plot of Fig.|2l 
the spectra decay linearly in the semi-log plot within the intermediate frequency regime, signi- 
fying that the modal correlation function is well-described by the algebraic function in the early 
stage. They gently decay at a high frequency, implying the existence of an extremely short-time 
dynamics (i.e., the bare-friction stage) dominated by the Markov part of the non-thermal fluctu- 
ations. Note that the power spectra of the KS equation exponentially decay in large-frequency 
regime 1I28] because the bare-friction stage does not exist there. 

Finally, the Fourier transform of the memory function due to the turbulent fluctuations is 



shown in Fig. [3l where Tk^u; was calculated by 

) 

drrfc(r) coswr. (23) 



The power spectrum r^^aj analytically relates to Uk,uj- Corresponding to Eq. (fTTI) . the evolution 
equation of Uk{t) can be represented as 

^ = - dt'r',{t - t'Wit') + R',{t), (24) 

where R'^it) denotes the non-thermal fluctuations and satisfies 

{R',it)ulito)) = 0. (25) 

Note that the integral range of Eq. (|24|) has been expanded compared to Eq. (fTTI) : nevertheless, 
the contribution is included in the fluctuating term BBIl . In the theoretical framework proposed in 
Ref. lai, the following relation holds; 

{R',it + r)R',{t)) 
^kiV = B ' ^^^^ 

indicating that F'^ ^ allows us to investigate the non-thermal fluctuations. It should be noted here 
that Eq. (|26l) holds in non-equilibrium systems. From Eq. (|24)) . the Fourier transform Uk,uj of Uk{t) 
satisfies 

(iw + Ffc [w] ) Mfc,,^ = i?^^^ , (27) 
where F^fcj] denotes the Fourier-Laplace transform of r'^(t), defined as 

POO 

r'fcM = / drrUr)e-^ (28) 



The power spectrum of Uk{t) is thus written as 



T[ 

\iu ■ 

Under the condition (fT2l) . Eq. (|29l ) reduces to 



Uk,. = ,. . .,2 - (29) 



= 7 ^\ , . (30) 

Although this is equivalent to Eq. (fT3l) . it might be more compatible with experiments such as the 
scattering methods. Because Tk,ui rapidly decay as shown in Fig.[3l Eq. (|30l) for a sufficiently large 
a; approximately reduces to 

10 



?7fc(r) = l-7f |r|, (31) 

i.e., the relaxation in the bare-friction stage is approximately linear decay, as mentioned in 
Ref. [29]. In contrast, at a sufficiently small to where Imrfc[t<j] — t- 0, we obtain the macroscopic 



friction coefficient [13 ih as 

(0) I (r) T / (u) 
Tk +% = l/^fc , (32) 

where rj^^ denotes the characteristic timescale of the modal autocorrelation function, defined as 

POO 

rf ) = / drUkir). (33) 
Jo 

IV. SUMMARY 

We have studied the relaxation dynamics of SMT, aiming to clarify statistical-physical prop- 
erties in the nonlinear nonequilibrium systems. Because the net autocorrelation function is a 
weighted mean of the modal autocorrelation functions, it has turned out that the CEF decay in 
SMT appears as a superposition of the dual structure. Further, this proceedings paper has concen- 
trated on the temporal power spectra of pattern dynamics and its modal element. The plot style of 
Fig. d is compatible with the viewpoint that the coexistence of the ballistic and diffusive motions 
produces a CEF decay. The physical origin of the relationship between CEF and complex systems 
such as the glassy dynamics is still an open question; however, the emergence of the CEF decay 
on SMT will provide us with a clue. In addition, we have shown the power spectra Uk,uj and Tk^^ 
and discussed the dual structure from a viewpoint of the harmonic analysis. Techniques such as 
the scattering experiments can measure 7*^°^ +rfc[a;] as the complex friction coefficients; therefore, 
the equations such as Eq. (|30l) will be helpful to investigate the turbulent fluctuations. 
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APPENDIX 



Soft-mode turbulence 



Instability in liquid crystal systems leads to pattern formation [|34l l35ll. Electrohydrodynamic 
convection is such a self-organizing phenomenon observed in nematic liquid crystal Oal . Let 
us consider systems of nematic liquid crystals with negative dielectric constant anisotropy ea = 
ey — e_i_, where ey denotes the dielectric constant parallel to the director n = 11(2;, y, z) representing 
the molecular orientation of a liquid crystal, e±_ the dielectric constant perpendicular to n. When 
the AC voltage larger than a threshold magnitude is a ppli ed to the systems, the electroconvection 



occurs due to the so-called Carr-Helfrich effect 1371 13811. interaction between the anisotropy of 
the liquid crystals and electrical current by impurity ions. The characteristic length scale of the 
electroconvection is much shorter than that of the heat convection, and consequently characteristic 
time scales of the electroconvection are easy to access in experiments. Thus, the electrohydro- 



dynamic convection has been studied as an ideal subject for the dissipative structure [|39ll . The 
squared magnitude V"^ and inverse frequency 1// of the AC voltage correspond to the Rayleigh 
number and the Prandtl number, respectively [40], where the former is proportional to the temper- 
ature difference and the latter is equivalent to the ratio of the viscous diffusivity to the thermal one 
in the Rayleigh-Benard convection. 

There are two types of layer alignment in nematic liquid crystal; one is planar alignment in 
which the director aligns parallel to substrates (x-direction), and the other is homeotropic align- 
ment in which the director aligns perpendicular to substrates (2;-direction). The rubbing to x- 
direction on the substrate's surface {x-y plane) is treated to make planer systems, and intrinsically 
breaks the rotational symmetry. In systems where nematic liquid crystals are homeotropically 
anchored at boundary {x-y plane), rotational symmetry are alive at the ground state as well as 
translation symmetry. With increasing magnitude V of the AC voltage parallel to 2;-axis, the 



Freedericksz transition occurs at a threshold voltage 114 111 . The directors n tilt in an arbitrary di- 
rection with respect to the 2-axis above Vy, accompanied by spontaneous violation of the rotational 
symmetry. The projection of the director onto the x-y plane, called C-director, can rotate on the 
x-y plane without requiring additional energy, i.e., C-director behaves as the Nambu-Goldstone 



mode 11421 - 14411 . With further increasing V , the electrohydrodynamic convection occurs at the elec- 



troconvective threshold voltage V^. The nonlinear coupling between the convective mode q(x, y) 

12 



and the Nambu-Goldstone one C(x, y) induces spatially and temporally disordered pattern called 



soft-mode turbulence (SMT) |19 



111) . SMT is regarded as an experimentally observed spatiotem- 



poral chaos. The diversity of the relaxation time at the onset of SMT implies that SMT arises 
supercritically at V^. The correlation length is finite but much longer than excited and dissipative 
lengths, where energy is injected at the excited length scale and dissipated at the smaller one. 

There are two types of SMT pattern; oblique rolls in / < /l and normal rolls in / > /l, where 
/l denotes the Lifshitz frequency l^llO]- The interaction between the convective and Goldstone 



modes are different between oblique and normal rolls [|45 



tend to become parallel to the C-director in the normal ro 



patch structure is observed only in the oblique rolls [113 



4611 ; the convective wave vector in SMT 
1 and oblique in the oblique roll. The 



We have studied the two-dimensio nal p attern dynamics of SMT. The experimental setup for this 



study refers to a standard one 



fly 



4811 ■ We filled the space between two parallel glass 



plates with the nematic liquid crystal, A^-(4-Methoxybenzilidene)-4-buthylaniline (MBBA). The 
surfaces of the plates were coated with transparent electrodes, indium tin oxide (ITO), which is a 
circular electrode with the radius 13 mm. The dielectric constant parallel to the director was 6.25 
and the electric conductivity parallel to the director 1.17 x 10^^ fi^^m^^. Note that the dielectric 
constant anisotropy is negative. The thickness between the plates is 50 fim for the net relaxation 
dynamics and 27 fim for the modal relaxation dynamics. To obtain the homeotropic alignment, 
the surface was laid by a surfactant, DMOAP (A^, A^-dimethyl-A^-octadecyl-3-aminopropyl- 
trimethoxysilyl chloride 50%). 

AC voltage V{t) = V2V cos{27t ft) was applied to the sample. The threshold voltage for 
electroconvection was 7.78 ± 0.05 V. We fixed / = 100 Hz that was much less than /l to obtain 
the oblique rolls. The temperature was set at 30.00 ± 0.05 °C and the magnetic field was absent in 
this study. Before each sampling, we waited for 10 min at Vw = 6.0 V and then for 10 min at set 
V, where Vp < Vw < K- The waiting time is sufficiently long for systems to become steady state. 

The transmitted light intensity at each pixel were digitized into 8-bit ( i.e., 256-level) informa- 
tion, where a series of pattern analyzing was processed according to Ref. l|49|] . 



Numerical algorithm to calculate memory function due to turbulent fluctuations 



In order to solve the memory function numerically by using an experimentally-obtained modal 
autocorrelation function U{t), we employed fitting by an expansion for short-time regime and 
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numerical integral of the evolution equation (U3l) for otherwise. Note that the wave-number de- 
pendence was omitted in the appendix; however, the coefficients actually depend on the wave 
number k. 

We first obtained the coefficients {a„} of an expansion around r = by fitting: 

oo 

f/(r) = l + ^a„r". (34) 

n=l 

Suppose that the memory function r(r) caused by the turbulent fluctuations is described by a 
continuous function, i.e., 

oo 

r(r) = ^6„r". (35) 

n=0 

Equation (fTSl) provides the bare friction as 

7^°) = -ai. (36) 
Moreover, the coefficients of the memory function as 

bo = ai^- 2a2 (37) 

and for n > 1, 

K = -{n + 2)(n + l)a„+2 + (n + l)aia„+i 
/ \ -1 

n 



m=l 



- ^ arahn-m | ) • (38) 

m 



For the long-time region, we solved the difference equation of Eq. (fTSl) . We first introduce 
A [/at as the central difference; 

Af/^ = ^^^^/^ ; ^^'^/^ (39) 
At 

where IJm = U (NAt), N denotes an integer, and At the time lag. Then, the evolution equation is 
digitalized to 

Un+1 — Un 



At 

-7^°^f>iv+i/2 - / dr'r{r')U{r-r'), 







(40) 
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where r should be interpreted as (A^ + 1/2) At. The convolutional integral is further digitalized to 

r dTr{r')U{T-r') 
Jo 



/ drr{r')U{T-T') 
Jo 

^ p{AI+l/2)At 

+ V / dr'r(r')f/(r-r') 

».r_i J(M-l/2)At 



At 



-7r'^oUN+l/2 + ^ A.tTMUN-M+l/2 (41) 



M=l 



with Fjv = r(A^At). Therefore, the memory function is solved numerically as 

2 UN+i-Ur, i^'Wn+i + Un 



N 



Ui + Uo A^ Ui + Uo 



2 f/i + f/o 



^ J2 Tm{Un~m+i + Un-m) (42) 



for > 1. The above discretization is not unique, of course; however, the way to discretize 
contributes only higher-order errors. 

The assumption (l35l) does not hold for r'(r) because it has singularity at r ~ originates from 
the Markov part of the non-thermal fluctuations. When one solves r'(r) numerically, the discrete 
analysis (|42l) with 7*^°^ = should be used for whole time region. Note that the initial value of the 
memory function is described by 



V - _ A^^^ (43) 
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